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I. INTRODUCTION 

The magnetization of type II superconductors is de- 
termined by the number of quantized magnetic vortices 
inside the sample [1] . As an externally imposed magnetic 
field is increased, vortices penetrate the sample in a pro- 
cess which at non-zero temperature is driven by thermal 
activation over energy barriers produced by the sample 
surface and by pinning centers in the bulk. When the 
external field is lowered, vortices leak out of the sample. 
The rate at which vortices move in and out of the sample 
determines the magnetic creep rate. 

Given that the magnetic relaxation is driven by ther- 
mal activation, it is rather surprising that experiments 
have found the creep rate to be essentially temperature 
independent in a wide temperature range [2-12]. Several 
mechanisms have been suggested to explain how com- 
peting factors are able to cancel the typical Arrhenius 
temperature dependence exp(C//T) for activation over an 
energy barrier U. The most prominent theoretical sug- 
gestion so far is probably the description in terms of col- 
lective creep [13]. 

Here we show that the lack of temperature dependence 
of the creep rate can be understood quite simply in terms 
of record dynamics[lA], which has recently been proposed 
as a general mechanism for the irreversible dynamics fol- 
lowing a sudden quench in glassy systems [14, 15]. 

We base our analysis on the following two experimental 
and numerical observations: 

A) For a long time glitches have been observed in 
the time dependences of the magnetic relaxation 
of type II superconductors (see [12] and references 
therein). As the external magnetic field is varied, 
the magnetization undergoes abrupt jumps when- 
ever vortices suddenly move in or out of the sample. 
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B) A number of studies, experimental as well as the- 
oretical [16-27], indicates that in a broad range of 
not too high temperatures the vortex system ex- 
hibits many characteristic features of glassy dy- 
namics. 

Many [28-35] but far from all [16-20] creep experi- 
ments study the response of the external magnetic field 
to a sweep, including sign reversal. For simplicity, we 
will presently concentrate on a setup in which the exter- 
nal field is initially ramped at a fixed rate up to a value 
and then remains constant for the entire duration of the 
experiment. 

We use a three dimensional version of the Restricted 
Occupancy Model (ROM) [21-27] to study the response 
to a fixed applied magnetic field. Important physical 
properties of the 3D layered ROM model have a nearly 
temperature independent dynamical evolution, matching 
in this respect simulation results previously obtained for 
the same model [21-27] as well as experimental results 
on type II superconductors [3-6, 11, 12, 16-20, 28-35]. 
As a strong temperature sensitivity is expected for the 
activated dynamics of an entirely classical model system, 
the mechanism behind the observed temperature inde- 
pendence warrants some further theoretical scrutiny. 

A fixed applied magnetic field can be expected to lead 
to a non-stationary time dependence of the internal mag- 
netization which slowly increases from zero to about the 
value of the applied field. Interestingly, the low temper- 
ature dynamical evolution of the response involves two 
types of configuration rearrangements, having widely dif- 
ferent time scales. One type consists of rapid glitches by 
which the magnetization irreversibly jumps to higher val- 
ues, as additional vortices enter the system. We refer to 
these glitches as quakes [15, 36, 37] to emphasize their 
non-equilibrium nature, their abruptness and their dra- 
matic effect on the state of the vortex system. The quakes 
are separated by much longer periods of apparent quies- 
cence, during which the vortex system is 'searching' for a 
configuration of larger stability and for a way to accom- 
modate more vortices inside the sample. Even though the 
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total magnetization docs not change, the internal spatial 
organization of the vortices undergoes considerable rear- 
rangements. 

Wc argue that the slow vortex creep associated with 
the quakes can be analyzed through the statistical prop- 
erties of record dynamics, which immediately explains 
the temperature independence of the creep rate[14] for 
the range of temperatures for which this description is 
applicable. By record dynamics wc mean the following. 
Consider a stochastic time signal xit)- The correspond- 
ing record signal is R{t) = max^<t{x(r)}, i.e. the largest 
value assumed by x(t) up to the current time t. The phys- 
ical idea behind record dynamics [14, 36, 38, 39] is that 
large irreversible configurational changes in noisy systems 
with a macroscopic number of metastable attractors are 
induced by noise fluctuations of record size. Similar be- 
havior is observed in other glassy metastable systems, 
e.g. gels and spin-glasses, and can be characterized sta- 
tistically in a similar way [15, 36, 37]. This highlights 
the underlying unity of non-equilibrium glassy dynam- 
ics at low temperatures and supports the possibility of a 
common theoretical description. 

The paper is organized as follows: The next Section 
summarizes the properties of the ROM model used in 
the simulations. Section III briefly introduces a possible 
mechanism [15, 38, 39] by which activated dynamics can 
become insensitive to the temperature in a glassy sys- 
tem, and demonstrates its relevance for the ROM model 
dynamics. Section IV focuses on creep rates and con- 
tains comparison with experiments. Finally, Section V 
presents a summary and a discussion. 



II. THE 3D ROM MODEL 

The two dimensional Random Occupancy Model 
(ROM) has been shown to reproduce the essential fea- 
tures of vortex dynamics at nonzero temperature [21 27]. 
Here we use Monte Carlo (MC) simulations of a general- 
ized three dimensional layered version of the ROM model 
to capture the long time relaxation of interacting vortex 
matter. 

In vortex matter, the length scale of the interactions 
can be very large compared with the average distance be- 
tween (pancake) vortices. At high densities, this implies 
that each vortex interacts with many others. For layered 
superconductors this situation can roughly be described 
by two length scales: the first is the range of the interac- 
tion parallel to the planes, this is the London penetration 
depth A. The second length scale is the vortex correla- 
tion length, ^||, parallel to the applied field (which we 
imagine to be perpendicular to the copper oxide planes 
for high temperature superconductors). The exact iden- 
tification of this length scale is difficult and is likely to 
depend on the anisotropy of the material, the nature of 
the pinning, the strength of the magnetic induction and 
on the temperature. This length scale may be related 
to vortex line cutting [40-4 7]. These length scales respec- 



tively give the horizontal and vertical lattice spacing of 
our model. The horizontal coarse-grained length scale 
corresponds to the penetration depth Iq = Xof the super- 
conducting material, and the spacing between the layers 
in our lattice we consider li ~ Smaller length scales 
are ignored. For our purposes this approximation is ac- 
ceptable because the length scales smaller than A seem 
to have little influence on the long time glassy properties 
of vortex matter. 

Another limitation of the model is that it ignores the 
variation of A with the temperature. As will become clear 
from our ensuing discussion of record dynamics, ignoring 
the temperature dependence of A is not crucial for our 
explanation of the observed temperature plateau of the 
creep rate. 

In a sample of a superconducting material the vor- 
tex matter behavior is determined by the competition of 
four energy scales [13]: intra and inter layer vortex- vortex 
interaction, vortex-pinning interaction and thermal fluc- 
tuations, all of which are schematically included in the 
ROM model. 

The Hamiltonian of the ROM model is thus the follow- 
ing: 

H = ^Aijninj-'^Aiini+'^A^ni+^ Ai (n, - rijf , 

ij i i {ij)^ 

a) 

where n; is the number of vortices on site i of the lattice. 
In a superconducting sample the number of vortex lines 
per unit area is restricted by the upper critical field (-Bc2) 
[1], so in the model the number of vortices per cell can 
only assume values smaller than Nc2 = Bc21q/4>o [25, 48]. 
Hence the name Restricted Occupancy Model. Moreover, 
as arc interested in a simiilation setup that does not 
require magnetic field inversion and the vortex-antivortex 
creation is strongly suppressed, we simply consider rii > 
0. 

The first two terms in Eq. (1) represent the repul- 
sion energy due to vortex- vortex interaction in the same 
layer, and the vortex self energy respectively. Since the 
potential that mediates this interaction decays exponen- 
tially at distances longer than our coarse-graining length 
A, interactions beyond nearest neighbors are neglected. 
We set Aii := Ao = 1, Aij := Ai if i and j are nearest 
neighbors on the same layer, and Aij := otherwise. 

The third term represents the interaction of the vor- 
tex pancakes with the pinning centers. A^ is a ran- 
dom potential and for the purposes of this work we con- 
sider that A^ has the following distribution P (Af) = 
{l-p)6 (Af ) - pS (Af - Ag). The pinning strength |Ag| 
represents the total action of the pinning centers located 
on a site. In the present work we use ]Aq] = 0.3. 

Finally the last term describes the interactions between 
the vortex sections in different layers. This term is a 
nearest neighbor quadratic interaction along the z axis, 
so that the number of vortices in neighboring cells along 
the z direction tends to be the same. 
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The parameters of the model are defined in units of Aq. 
The time is measured in units of full MC sweeps. The 
relationship between the model parameters and material 
parameters is discussed in [25, 48]. 

Each individual MC update involves the movement to 
a neighbor site of a single randomly selected vortex. The 
movement of the vortex is automatically accepted if the 
energy of the system decreases; if the energy of the sys- 
tem increases, the movement is accepted with probability 
exp(-A£;/T) [49]. 

Given that the movement of pancake vortices is re- 
stricted to the superconducting planes we only allow MC 
movements parallel to the planes. We have used periodic 
boundary conditions along the z direction. 

The external magnetic field is modeled by the edge 
sites on each of the planes. The density at the edge is kept 
at a controlled value. During a MC sweep vortices may 
move between the bulk sites and the edge sites. After 
each MC sweep the density on the edge sites is brought to 
the desired value. Initially the external field is increased 
to a desired value {Ng^t = 10 vortices per edge site) by a 
very rapid increase in the density on the edge sites. We 
have here used a sweeping rate 7 of 0.25 per MC sweep 
(compared with 7 G [10~^, 10^^] in our previous studies 
[21-27]). After this fast initial ramping the external field 
is kept constant, while we study how the vortices move 
into the sample. The age of the system, tw, is taken to 
be the time since the initial ramping. 

We have studied systems of different sizes, and obtain 
similar results except for very small system sizes. Our 
key results were obtained in a system consisting of 8 lay- 
ers of size 16 X 16. The model parameters used in our 
simulations were: 

number of realization = 4000 
Ai = 0.28 
A2 = 0.5 
p = 0.5 
AP = 0.3 
Nc2 = 27 
7 = 0.25 

A^ext = 10 

III. RECORD DYNAMICS 

We will in this section describe how the observed 
plateau in the temperature dependence of the magnetic 
creep rate can be explained in the framework of record 
dynamics. Let us first mention the salient features of 
what we mean by record dynamics. Consider a stochas- 
tic signal x{t) with no time correlations. Now derive 
the record signal R{t) = maxT-<t{x(r)}. We note that 
R{t) is a monotonous piecewise constant function which 
only increases its value at discrete times tk, whenever 
x{t) manages to fluctuate to a value larger than any en- 
countered previously. For our present purpose, the most 



important property of the statistics of the record times 
tk is that the probability that exactly q records occur 
during the time interval [tw, tw + 1] (where is the time 
since the initiation), is to a good approximation Poisson 
distributed on a logarithmic time scale [14, 36], i.e., 

p(g) = M!exp(-(g)) (2) 

with average number of quake events proportional to log- 
arithmic time 

(g) =alog(l+i/i^). (3) 

Here a is the logarithmic rate of events. To get the gist 
of the mathematics behind Eq. (2) (for full detail see 
ref. [14]), we note that since the largest outcome, i.e. the 
record, of t independent trials is equally likely to occur 
at any of the t instances in the sequence, it occurs at the 
first attempt with probability 1/t. Hence, the probability 
p(l) of exactly one record in t trials is 1/t. independently 
of the distribution of the underlying signal x{t)- It is im- 
portant to point out that the same also holds for the gen- 
eral expression for p{q) in Eq. (2). The independence of 
p{q) on the distribution of random numbers corresponds 
to the independence of x(t) on the thermal noise and will 
translate directly into the temperature independence of 
the creep rate. 

It has recently been shown [15, 38, 39] that in (glassy) 
systems having a large number of dynamically inequiv- 
alent attractors, temperature independence of suitably 
coarse-grained dynamical variables can arise from the pe- 
culiar way in which the attractors are selected as the sys- 
tem evolves from a typically rather unstable initial con- 
figuration through gradually more stable ones. A similar 
noise insensitivity of stochastic dynamics has been ob- 
served with other types of noise, e.g. in driven dynamical 
systems [14, 50] and in evolutionary dynamics[15, 51]. 

That the ROM model also exemplifies this type of be- 
havior can be gleaned from Fig. 1, showing the time de- 
pendence of the number of vortices N{t) which during 
a single run have entered the system up to time t. Im- 
portantly, the length of the quiescent periods typically 
increases with time - notice the logarithmic time axis 
in Fig. 1. Were this not the case, the dynamics would 
appear continuous in terms of a suitably coarse-grained 
time scale. Conversely, the lengthening of the intervals 
between successive quakes signals the anticipated grad- 
ual entrenching of the dynamics into dynamically more 
stable configurations. 

Also important is that the overwhelming majority of 
the observed glitches lead to states with a higher num- 
ber of vortices. This de-facto irreversibility of the dy- 
namics enables one to meaningfully approximate the sig- 
nal N{t) with the monotonically increasing record sig- 
nal R{t) = max[o,f]{A^(t)}. We stress that this theoreti- 
cally convenient idealization is only applicable within the 
strongly non-equilibrium regime of our present concern. 

Statistical insight into the time evolution in the num- 
ber of vortices present within the system is provided by 
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FIG. 1: The detailed time variation of the total number of 
vortices N(t) on the system for a single realization of the 
pinning potential and the thermal noise in a 8 x 8 x 8 lattice 
for T = 0.1. Notice the monotonous step function character 
of the time series. 




N(t+t )-N(t ) 

FIG. 2: (Color online) The main panel contains the temporal 
evolution of the probability density function, P{N{t + t^) — 
N{tw)), of the number of vortices entering for ~ 1000 and 
three different observation times t = 188, 2791, 8371 given by 
the black, slightly jagged, curves. The smooth curves (blue) 
are a fit to the theoretical expression (see Eq. (4)). The 
system is 16 x 16 x 8 and T = 0.0001. The insert shows 
the quake size distribution for various temperatures for the 
time interval between t = 1000 and t = 10000. For T < 0.1 
the distribution has an approximately exponential tail. For 
T — 0.5 the shape gets closer to a Gaussian. 



Fig. 2, where the empirical distributions of N{t) are dis- 
played for three different times, which are equidistantly 
placed on a logarithmic time axis. The insert in Fig. 2 
shows the tail of the probability density function (pdf) 
of the number of vortices, p{v), entering during a single 
quake. To a good approximation, the tail is exponential 
and the observed time dependence and temperature de- 
pendence of of p{v) are negligible, except for the highest 
temperature T — 0.5. 

The interpretation in terms of record dynamics sug- 



gests that the probability that exactly q quakes occur 
during the time interval [1^ + t\ is Poisson distributed 
on a logarithmic time scale [14, 36] according to Eq. (2). 
We approximate the pdf for the number of vortices v 
which enter during a given quake event (see insert Fig. 
2) by an exponential distribution p{v) — exp{—v/v)/v, 
and assume that subsequent quakes are statistically in- 
dependent. The number of vortices entering during ex- 
actly q quakes is then a sum of exponentially distributed 
independent variables, and is hence Gamma distributed. 
We finally obtain the pdf for the total number of vortices 
entering during [tw +t] by averaging the Gamma dis- 
tribution for q quakes over the probability Eq. (2) that 
precisely q quakes occur within the time interval of inter- 
est. This leads to the following expression for the pdf of 
total number of vortices AiV = N{t -f iu,) — N{tw) which 
may have entered during the time interval [tm ,tw + t] 



p{AN,t)^ 




where /i denotes the modified Bessel function of order 
1[52]. The above theoretical prediction, is compared in 
Fig. 2 with our simulation results. To estimate (q) ac- 
cording to Eq. (3), we used a = 22.6, as obtained from 
the logarithmic rate of the quake events. We can deter- 
mine the average v in two ways. Either directly from the 
simulated distributions in the insert of Fig. 2 or from fit- 
ting Eq. (4) to the simulated data in the main frame of 
Fig. 2. In both cases we find v — 16. We also find that v 
is essentially temperature independent for temperatures 
below T « 0.1. This is also expected from the insert 
in Fig. 2. It is important to mention that the MC dy- 
namics does overcome plenty of positive energy barriers, 
AE > 0, through thermal activation for temperatures 
in the range 0.01 < T < 0.1. As the temperature is 
lowered fewer MC updates correspond to AE > and 
for the lowest temperatures MC steps involve AE < 
only [24]. Nevertheless, the record dynamics remains es- 
sentially temperature independent for T < 0.1. 

The agreement is encouraging and suggests that the 
process of vortex penetration into the sample can be de- 
scribed in terms of a Poisson process with logarithmic 
time argument, for short the log-Poisson process. We 
also note that the log-Poisson statistics covers the tem- 
poral distribution of the quakes but has nothing to say 
on the size distribution of the jumps, i.e. the number v 
of vortices entering during a single event. This stochastic 
quantity could in principle introduce additional time and 
temperature dependencies. However, as mentioned, the 
insert of Fig. 2 shows that for a very broad parameter 
range this is not the case. Accordingly, the creep rate ob- 
tained by convoluting the distribution of the quake sizes 
and the log-Poisson distribution of number of quakes will 
also be temperature independent. 

The link between record statistics and the stochas- 
tic dynamics of a glassy system is discussed in detail 
in ref.[36] on the basis of several idealized physical as- 
sumptions. The first element is the existence of a large 



5 



number, in principle a continuum, of attractors. These 
are sets of configurations clustered around a local en- 
ergy minimum and supporting equilibrium-like reversible 
thermal fluctuations. By contrast, attractor changes — 
our quakes — are assumed to be irreversible on the time 
scale at which they occur. The exact nature of the quakes 
is not entirely clear in our system. At intermediate tem- 
peratures they are related to activation over barriers and 
the jump in N{t) is associated with a increase in the en- 
ergy of the interacting vortices. However, at the lowest 
temperatures there is not sufficient thermal energy avail- 
able for the system to climb over energy barriers. The 
thermal fluctuations are only able to push the vortices 
along equipotential trajectories or to lower potential en- 
ergy configurations. In this regime the vortex motion is 
hindered by jamming and the quakes are of a mechanical 
nature [24]. 

An interpretation in terms of record dynamics implies 
that the dynamical bottlenecks overcome by fluctuations 
are determined by the actual noise history, and not pre- 
determined in a static fashion. For other model sys- 
tems [14, 50], the validity of the assumed linkage between 
noise records and barriers was confirmed by considering 
white noise perturbations drawn from a distribution with 
finite support, e.g. a box distribution, and by then study- 
ing the properties of the selected attractors as a function 
of the maximum size of noise. 

Record-induced dynamics has thus a number of 
testable predictions, the most interesting of which is, 
for our purposes, the logarithmic time dependence and 
the striking temperature independence of the number of 
quakes occurring in the time interval [0, i], which are 
shown in the following Section. 

IV. CREEP RATES 



Let us now turn to the dependence on time and tem- 
perature of the total number of vortices in the sample. 
At time t — we rapidly increases the external field from 
zero up the value iVea;t (see section II). The vortex den- 
sity of the bulk sites gradually increases as vortices move 
in from the boundary. 

In Figure 3 we present, for a very broad range of tem- 
peratures, the average density n{t) of the bulk sites as 
function of the natural logarithm of time log(<:). As antic- 
ipated, the time dependence is temperature independent 
for all but the highest temperatures. 

One can identify three different temporal regimes sep- 
arated at times ti « 300 and t2 ~ 3 x 10*. For the re- 
maining of this paper we will focus our analysis on the in- 
termediate regime ti < t < t2 (and choose tw = ti). Our 
reason for this is that at ti vortex interactions become 
essential through out the entire system. The late time 
regime t > t2 is very difficult to resolve appropriately in 
simulations and probably equally difficult to study ex- 
perimentally. 
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FIG. 3: (Color online) Vortex density time dependence in 
the system for various temperatures. For T < 0.1 the vortex 
number is a piecewise linear function of log(f). For T = 0.5 
the system relaxation becomes faster. In this plot n{t) is 
averaged over many system realizations. 

For times ti < t < t2 Fig. 3 demonstrates that N{t) 
depends linearly on log{t) to a very good approximation. 

The linear logarithmic time dependence is of course 
entirely consistent with the record dynamics outlined in 
the previous section. We consider the total number of 
vortices in the system N{tw + 1) to be the accumulated 
effect of vortices entering during quake events that have 
occurred prior to time + t. Let denote the time 
of occurrence of quake number k and let vu denote the 
actual number of vortices entering during this quake. We 
then have 

N(t + t^)^N{t^)+ Vk, (5) 

t^u<tk<tu,+t 

where the sum is over all quakes that occured during the 
time interval [tw ,tw + t] 

From Fig. 2 we know that Vk is temperature inde- 
pendent and possesses a well defined average v. Since 
the average number of quakes increases according to Eq. 
(3) , record dynamics predicts the following (temperature 
independent) temporal evolution of the average number 
of vortices 

AiV = {N{t + t^)) - {N{ty,)) - av log(l + t/ty,). (6) 

i.e. for the considered time regime t/tw ^ 1 a tempera- 
ture independent logarithmic rate given by 

dAN/d log(t) ^ av/{l + tn,/t) ^ av. (7) 

We extract the rate of the quakes in the simulations 
from temporal signals like the one exhibited in Fig. 1. 
In Fig. 4 we demonstrate that the quake rate is indeed 
approximately independent of temperature in the broad 
temperature interval lO^'* < T < 2 x 10~^. 

Finally, Fig. 5 shows the near temperature indepen- 
dence of the actual rate with which N{t) changes. We 
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FIG. 4: Quake dynamics rate, a, versus T for the time interval 
between t — 1000 and t = 10000. Notice that a is roughly 
temperature independent for three orders of magnitude below 
T < 0.1. 
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FIG. 5: (Color online) Main panel: Numerical results for the 
creep rate versus T for the time interval between t — 1000 and 
t = 10000. In the low temperature region the creep rate is 
constant within our numerical precision for about two orders 
of magnitude - we observe a nonzero creep rate in the T 
limit. Insets: experimental results for the creep rate versus 
T. The right inset shows data from Keller et al. for melt pro- 
cessed YBCO crystals with the magnetic field applied along 
the c axis (squares) and ab plane (circles). The left inset shows 
data from [53] for unirradiated (squares) and 3 MeV proton- 
irradiated (circles) YBCO flux grown [54] crystals with a 1 T 
magnetic Field applied parallel to the c axis. 

extract this rate from the data shown in Fig. 3 in 
the time region ti < t < t2 and plot the normalized 
creep rate S = d\og[M{t)]/d\og{t) in order to com- 
pare consistently with experiments, here we have used 
Mit) = \Nit)-N,,t\. 

V. CONCLUSION 

We have presented an analysis of simulated vortex 
creep data in terms of record dynamics. This approach 



allows us to interpret the observed temperature indepen- 
dence of the creep rate as a generic property of the dy- 
namics of records obtained from the underlying fluctuat- 
ing sequence. To establish the temperature independence 
of the creep rate we do not need to know the detailed na- 
ture of the quantity being gradually maximized. Nor do 
we need a description of the intermittent vortex quakes 
which are responsible for the abrupt changes in the num- 
ber of vortices. All we make use of is the assumption, sup- 
ported by the simulated model, that the abrupt glitches 
in the number of vortices inside the sample can be in- 
terpreted as arising from the records of some stochastic 
process. We showed that the simulated creep rate be- 
haves in a way very similar to published experimental 
data for YBCO. 

It is obvious that a better experimental and theoreti- 
cal understanding of the nature of the vortex quakes is 
interesting and future study of the ROM model will seek 
to improve our understanding of the spatial and dynam- 
ical properties of these quakes. We can already conclude 
that the physical mechanisms involved must be different 
at low and high temperatures. At the lowest temper- 
atures activation over free energy barriers are excluded 
and the quakes are related to mechanical rearrangements 
of the vortices[24]. At elevated temperatures the quakes 
are expected to be triggered by activation over thermal 
barriers. The present paper shows that record dynam- 
ics can be used to understand the temperature range 
from very low temperatures, where no barriers can be 
climbed, up through a regime where thermal activation 
does take place. For high temperatures (in our case for 
T > 0.1) the description in terms of record dynamics 
breaks down. This happens when there is sufficient ther- 
mal energy available to make any trapped metastable 
configurations short lived. 

Let us finally mention that our description in terms of 
record dynamics may not exclude aspects of previous de- 
scriptions of vortex relaxation in terms of e.g. correlated 
collective vortex creep [13]. We would rather think of 
our approach as contributing to an understanding of the 
detailed nature of the dynamics of the correlated vortex 
regions. 
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